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£-( Abstract 

The process of neutrino-pair radiation due to formation and breaking of 



^NJ ■ Cooper pairs of protons in superconducting cores of neutron stars is considered 

>: 

with taking into account of the electromagnetic coupling of protons to ambient 
electrons. It is shown that plasma polarization strongly modifies the effective 

o' 

^^ ■ vector weak current of protons. Collective response of ambient electrons to the 

O 

c| ■ proton quantum transition contributes coherently to the complete interaction 

a: 

with the neutrino field and enhances the rate of neutrino-pair production by 



-C ■ two orders of magnitude. 
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Except for the early stages of evolution, temperatures inside neutron stars are lower 
than the critical temperatures for proton superconductivity T cp and/or neutron superfluidity 
T cn . Conventionally, neutron superfluidity as well as proton superconductivity, exponentially 
reduce the heat capacity of neutron matter and the rates of neutrino production from nucleon 
reactions in the core of a neutron star |TJ, , 0, ||. When T < T cp ,T cn , nuclear matter 
in the stellar core consists of two components. One part of nucleons forms a condensate of 
Cooper pairs, the other one represents unpaired quasi-particles. Both processes, Cooper- 
pair formation and pair-breaking coexist in statistical equilibrium and result in additional 
neutrino-pair emission from the neutron star. This mechanism of energy loss from superfluid 
neutron stars was first proposed and evaluated for neutron singlet-state superfluidity by 
Flowers et al. |J, but it has been included into cooling simulations only recently |7|], || , |§. 
It was shown, that the phase transition of neutrons in a superfluid state does not decelerate 
the process of cooling, but, on the contrary, neutrino emission due to formation and breaking 
of Cooper pairs accelerates the neutron star cooling. 

Cooper pairing of protons take place likely in 1 Sq state |J. Since the total spin of the 
Cooper pair in the 1 Sq state is zero, the axial- vector contribution of the weak neutral current 
vanishes. Contribution to vv emissivity comes only from the vector weak interaction. There- 
fore, neutrino emission produced by proton pairing is conventionally assumed to be strongly 
suppressed due to the numerical smallness of the weak vector coupling constant of protons 
|T0| . Such inference is made on the basis of calculations which ignored electromagnetic cor- 
relations among the charged particles in a QED plasma. Actually, protons in the plasma 
are coupled to ambient electrons via the electromagnetic field. By undergoing a quantum 
transition to the paired state, protons polarize the medium, thus inducing the motion of 
electrons inside the Debye sphere around them. The electron weak current associated to 
this motion generates neutrinos coherently with the weak current of protons, because the 
wavelength A of radiated neutrino pairs is much larger than the electron Debye screening 
distance D e (typically, D"l/\ 2 < 1CT 2 ). In the present article we study an effective weak 
current of protons with taking into account the collective contribution of ambient electrons, 



and estimate the neutrino-pair emissivity due to formation and breaking of proton Cooper 
pairs. 

In what follows we use the system of units h — c — 1 and the Boltzmann constant 
ks — 1- The low-energy Lagrangian of proton interaction with the neutrino field in vacuum 
can be written in a point-like approach 

G F _ 

Cvac = 7K U 7n (1 - 7s) »f, (1) 

where Gp is the Fermi coupling constant, and the proton weak current is of the form 

f = ^{Cy - C Al ,)^ (2) 

Here ifj stands for the proton field; Cy and Ca are the vector and axial-vector coupling 
constants of the proton, respectively. Reserving the capital letter notations Cy and Ca 
for proton coupling constants, we will use, at the same time, small letters, Cy, for electron 
coupling constants with cy = | + 2sin 2 9\y ~ 0.96 for electron neutrinos, and c' v = —\ + 
2 sin 2 9w — —0.04 for muon and tau neutrinos; 9w is the Weinberg angle. 

To describe the above-mentioned collective effects, one has to include in the matrix 
element of the reaction two Feynman diagrams shown in Fig. 1. The sum of these diagrams 
represents the effective weak coupling of in-medium protons with the neutrino field. The 
matrix element of this interaction should be calculated between exact initial and final states 
of two proton quasi-particles undergoing quantum transition to the paired state. The first 
diagram in Fig. 1 represents the proton weak interaction in vacuum, and the second one is 
the contribution of ambient electrons. Some explanations are necessary here. At the first 
sight, the second diagram contains additionally a small factor e 2 = 1/137. However, the 
fine structure constant enters the matrix element in a combination with the electron number 
density, known as the Debye screening distance. Therefore, the parameter of the problem 
is k 2 D 2 , where k is the momentum carried out by the neutrino-pair. As will be shown, 
in the case k 2 D 2 < 1 the second diagram gives an important contribution and should be 
necessarily included in the matrix element of the reaction. 



The electron loop in the second diagram can be expressed in terms of polarization tensors 
of the electron gas, defined in the one-loop approximation 

U^(K)=4me 2 TrJ-0^YG(p)YG(p + K), (3) 

nf (K) = 4vne 2 Tr J ^- Y G(p) Y l5 G(p + K). (4) 

Here G(p) is the in-medium electron Green's function, and K = (u, k) is the four-momentum 
transfer. Thus, the in-medium effective weak current of protons is of the following form 

JM = ^{C v - CU 75 )V - -^ (^tV) D Xp (evil'" - c A Iit) , (5) 

where D\ p (K) is the photon propagator in the medium. An extra factor (Air) appears 
here because the factor 4n is traditionally included in the definition of polarization tensors. 
The minus sign in front of this factor is due to the fact that the second diagram includes 
a proton charge e > in the left-hand electromagnetic vertex of the virtual photon, and 
an electron charge — e in the right-hand vertex. The corresponding e 2 factor which appears 
accordingly to the second diagram, is also included in the definition of polarization tensors. 
To write the in-medium photon propagator one has to take into account that the super- 
conducting electric current of protons and the normal electric current of electrons coexist, 
participating in electromagnetic oscillations of the medium. The total Lagrangian density 
of the complex field \I/ of proton Cooper pairs, and electrons interacting with the electro- 
magnetic field Afj_ has the form 

L=\(d, + 2teA,)y\ 2 + M 2 Cp \y\ 2 -\\y\ 4 --^F* u -j;A» + L° e (6) 

where A is a constant for proton-proton interaction resulting in the proton pairing, and 
Mc p — 2M* is the mass of the Cooper pair consisting of two protons of effective mass M*; 
Fpv = d ii A v — d u A^ is the tensor of electromagnetic field; j£ is the electron current, and 
L° is the Lagrangian density of free electrons. As a result of spontaneous breaking of the 
symmetry of the ground state of the system, the vacuum expectation value of the Cooper 



pair field is nonzero ( |\l/| ) =^l- Since the ground state of the system corresponds to zero 
spin and zero total momentum of the Cooper pair the vacuum expectation value is connected 
with the number density of Cooper pairs iV s by the relation N s = 2E^l, where E ~ 2M* 
is the energy of the Cooper pair corresponding to zero total momentum, and N s is one-half 
the number density of paired protons N p . Thus, 

^o = — ■ (7) 

8M; v ; 

We represent the field of the Cooper condensate in the form 

$ = $ o (l + p)exp(-2ie0) (8) 

where p and <fr are arbitrary real functions of the coordinates and time, and subsitute in the 
Lagrangian density (|G]) for the fields p and <fi. Taking into account the following identity 

\{d„ + 2ieA ll ) ^| 2 = ^l (d,p) 2 + 4e 2 ^ 2 (A, - d^f (1 + pf (9) 

we can make the gauge transformation A'^ = A^ — d^cj). Since the quantity F^ u = F' = 
d^A' v — d u A' as well as the electron current j e = j' e are gauge-invariant we obtain 

L = *P 2 o {d,pf + Ae^lA'l (1 + pf + M 2 p |^| 2 - A |^| 4 - ^F% - ftA'" + L\ (10) 

As a result, the Goldstone field <fi is absorbed by the gauge transformation. Considering 
this and taking variations of the Lagrangian with respect to the field A' 11 , in the linear 
approximation we obtain the equation 

dJT^ + 327^^ = 47TJ-, 0^ = (11) 

In the absence of the electron current (j* = 0J, this equation would describe the eigen 



photon modes of mass m 7 = W327re 2 \l/ 2 1 = J4:irN p e 2 /M* This is the well-known Higgs 
effect. However, the induced electron current should also be taken into account. In the 
momentum representation Aiij e = —U^A'", and the Eq.(|ll|) takes the form 

^;-mX-n F (^)^ = o, i^4 = o (12) 



To specify the components of the polarization tensor, we select a basis constructed from 
the following orthogonal four-vectors h p = (u;,k) /\fK 2 , l^ = (k,un) /VK 2 , where 

the space-like unit vector n = k/fc , k — |k| is directed along the electromagnetic wave 
vector k. Thus, the longitudinal basis tensor can be chosen as L pp = —l p l p . The transverse 
(with respect to k) components of Tl pp have a tensor structure proportional to the tensor 
T Pfl = (g pp — h 9 ^ + l p l p ) , where g ppl = diag(l, — 1, —1, —1) is the signature tensor. In this 
basis, the polarization tensor has the following form 

Y[p p (K) = in(K)L pp + n t (K)T pp , (13) 

In the case of a strongly degenerate ultrarelativistic electron plasma (vp ~ 1), the longitu- 
dinal and transverse polarization functions areQ ||11|| : 

-M-OHM 1 ^))- (15) 

The electron plasma frequency and the Debye screening distance are defined as 

2 4 2 2 47rn e e -^ o 2 rta\ 

Se = Yn 6 "■ = -^-' Df = ^ (16) 

with /i e and n e being, respectively, the chemical potential and the number density of elec- 
trons. The axial polarization tensor is given by 

nr (K) = n A (K) ih x eP^°, n A (K) = ^—^ (K) . (17) 

The total energy of neutrino-pairs is uj ~ A <C /i e . When u > k, polarization functions 
differ in the order of magnitude, namely, tta ~ 7Tj,tA/jU e . Therefore, the axial polarization 
of the medium can be neglected. 



1 Our Eq. ( fUD differs from Eq. (A39) of the Ref. [11] by an extra factor (u> 2 /k 2 — l) because 
our basis l p , h p is different from that used by Braaten and Segel. All components of the complete 



tensor ( |l3|) identically coincide with that obtained in [11] for the degenerate ultrarelativistic case 



In the absence of the electron current, the propagator of a massive photon is of the form 

4n ( K^K V \ 

D o v = 1 7f L -lU v -^-) (18) 

The in-medium photon propagator, taking into account the induced electron current, can 
be found with the aid of the Dyson's equation. The solution is: 

471 A7T 4:71 

^w- p-kV " ^-^-/ "'^ (19) 

Due to conservation and gauge invariance of the electron current hxll^ = Yi x ^h^ = 0, the 
last term of Eq.( [19|) does not contribute to the effective weak current of a proton in the 
medium. 

Insertion of (pT3|) and (pT9|) into (|5|) with the above approximation yields the effective weak 



current of a proton in the medium. 

J" = ^f{C V - C Al5 )^ - Cy (farf) ( 2 J'^, + 2 ^f X \ ) (20) 

v ' \ oo z — fir — m z — txi ui z — k z — m z — n t J 

The poles of this expression give the eigen photon modes in the medium. In the case of 
K 2 > the eigen mode frequency u (k) is larger than the plasma frequency of electrons. 
At temperatures T < T c <^. u pe , the number of such excited oscillations in the medium is 
exponentially suppressed, therefore contributions from the poles can be neglected. 

The contributions of longitudinal and transverse virtual photons to the effective weak 
current (|2CiD are proportional to the following factors 

Fl (k : u)= ^ ^ (21) 

to 2 — k 2 — m 2 — iii D 2 uj 2 (1 — x 2 ) — m 2 D 2 + cp(x) 

F (k n)= ^ = l + t ( x ) t^\ 

U ' } ~ u 2 -k 2 -m 2 -7i t " 2D 2 e u 2 {l~x 2 )-m 2 D 2 -(l + ip(x)) l ' 

with x = k/u and 

1 i„ 1 + x \ 2 n 2 _ Me 



^^-.(l-^l-ih.lil , «W = £; (23) 



The latter equality takes into account local neutrality N p = n e of the medium^. 

The neutrino wave-length is much larger than the Debye screening distance, therefore, 
Dluj 2 (1 — x 2 ) < 4A 2 _D 2 (1 — x 2 ) <C if (x). Neglecting this small term we obtain 

*-,(») = - ,f'*' ■ F t{x) = - —l±^l— 2 . (24) 

Lp (x) - mi.D 2 e 1 + Lp (x) + mtp% 

Due to the electron contribution the effective proton interaction with the neutrino field has a 
complicated form. According to Eq.(|20|), the vector part of the effective proton weak current 
is substantially modified by polarization of the medium. Assuming the proton effective mass 
to be M* ~ 0.8M P , functions Fi (x), F t (x) are plotted in the Fig. 2 for typical values of the 
electron chemical potential /i e ~ 100 MeV and fi e ~ 200 MeV. 

For the moment, we shall estimate the contribution of the medium polarization to the 
neutrino emissivity caused by formation and breaking of proton Cooper pairs and leave the 
complete calculation for a future work. To roughly estimate the contribution of collective 
effects, we take advantage of the fact that Fi(x) and F t (x) are slowly varying functions, and 
replace them in Eq.(|20"D by a constant value —0.8. With this simplification we obtain 



J» ~ ^[C v - C A 1 5 )^ + 0.8c v (^7a^) (L x » + T A ^) . (25) 

Taking into account the identity L Xfl + T Am = g x ^ — /i A /i M , and gauge invariance of the 



electromagnetic current of protons (ip , y\ip) K = 0, we can reduce the effective proton 
weak current ( pq) as follows 

r ~ i)^{C v + 0.8c y - C A75 )^. (26) 

Thus, the effective vector weak coupling of in-medium protons is of the order of Cy + 0.8cy , 
in contrast to Cy = ~ — 2 sin 2 9\y ~ 0.04 for bare protons in vacuum. If protons interact with 
electron neutrinos then Cy + 0.8cy ~ 0. 8 , while for muon and tau neutrinos Cy + 0.8c' v ~ 
.008. The total spin of a Cooper pair in the 1 5'o state is zero therefore, the contribution to 



Strictly speaking, this equality is valid only at zero temperature, when all protons are paired. 



vv emissivity comes only from the vector weak interaction. In previous calculations ||, |10 
the vv emissivity was proportional to C v . Taking into account the collective effects, one 
has to replace Cy by Cy + 0.8cy for electron neutrinos and by Cy + 0.8c' v for muon and tau 
neutrinos. Thus, the vv emissivity is enhanced about two orders of magnitude: 



Q (C v + 0.8c v y + 2(C V + 0.8c' v 
Qo 3C v 



— — ~ 133. (27) 



This ratio demonstrates the large importance of collective processes in the core of a super- 
conducting neutron star. As Cy <C cy, the leading contribution to Q comes from the weak 
current of electrons which is electromagnetically induced by the quantum transition of the 
initial proton. 

The neutrino emissivity due to Cooper pairing of protons, as well as that for neutron 
pairing is proportional to the Fermi momentum of the relevant particles, i.e. it only weakly 
( ~ nl/n) depends explicitly on the partial number density. Besides that, both the superfluid 
and the superconducting energy gaps, calculated theoretically, sensitively depend on the 
model of strong interactions and vary in the range from ten keV to some MeV JIBf , | |14| , 



[I5fl, flTTf]. Thus, though the number density of protons is much smaller than that of neutrons, 



their partial contribution to the total energy losses is very sensitive to the ratio between the 
energy gaps of neutrons and protons. It can be dominant at some cooling stage, depending 
on temperature, T cn , T cp , and density. Therefore, the neutrino radiation, caused by a 
singlet-pairing of protons, should be taken into consideration in the scenario of cooling 
of superconducting neutron stars at the same level, as other basic processes of neutrino 
production. 

We have shown in the present paper that the effective vector weak current of in-medium 
protons is strongly modified by ambient electrons. Their weak current, electromagnetically 
induced by the proton quantum transition, enhances neutrino production due to Cooper 
pairing of protons by two orders of magnitude. An analogous conclusion can be made about 
neutrino-pair production due to Cooper pairing of H - hyperons. Though the vector constant 
of the weak coupling for S~ hyperons is small relatively to those for neutrons and S ± , the 



interaction of S with neutrinos is enhanced in the medium by ambient electrons. 
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Fig 1. Feynman graphs contributing to the effective proton interaction with a neutrino 
field in the medium. The first diagram presents the vacuum-like weak interaction. The 
second diagram, with the electron loop, describes the interaction via an intermediate virtual 
photon, shown by the dashed line. 
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Fig 2. Functions Fi (x), F t (x) are plotted for two cases of the electron chemical potential 
\i e ~ 100 MeV (deshed lines) and /i e ~ 200 MeV (solid lines). The proton effective mass 



M* = 0.8M P 
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